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Abstract — -Waveguide Meshes are efficient and versatile 
models of wave propagation along a multidimensional ideal 
medium. The choice of the mesh geometry affects both 
the computational cost and the accuracy of simulations. In 
this paper, we focus on 2D geometries and use multidimen- 
sional sampling theory to compare the square, triangular, 
and hexagonal meshes in terms of sampling efficiency and 
dispersion error under conditions of critical sampling. The 
analysis shows that the triangular geometry exhibits the 
most desirable tradeoff between accuracy and computational 
cost. 

Keywords — Waveguide meshes, finite difference methods, 
wave propagation, dispersion error, multidimensional sam- 
pling. 



I. Introduction 

AMONG the techniques for modeling wave propaga- 
tion in multidimensional media, the Digital Waveg- 
uide Meshes have recently been established 0] , |Q , ||] , Q , 
Hi Hi 7 @ as intuitive and efficient formulations of finite 
difference methods ||. 

A Waveguide Mesh (WM) is a discrete-time computa- 
tional structure that is constructed by tiling a multidimen- 
sional medium into regular elements, each giving a local 
description of wave propagation phenomena. This local de- 
scription is lumped into a waveguide junction ||, which 
is lossless by construction. Therefore, waveguide meshes 
are free of numerical losses, even though lumped passive 
elements can be explicitely inserted to simulate physical 
losses. However, wavefronts propagating along a multidi- 
mensional WM are affected by dispersion error, i.e. differ- 
ent frequencies experience different propagation velocities. 
Numerical dispersion cannot be completely eliminated, but 
it can be arbitrarily reduced increasing the density of the 
elements, and minimized choosing the "least dispersive ge- 
ometry". Moreover, interpolation schemes or off-line 
warping techniques JHJ can be applied to attenuate the 
effects of numerical dispersion. 

In this paper we investigate how the density of waveg- 
uide junctions and the sampling frequency affect the signal 
coming out from the model. As the analysis depends on 
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the geometry of the WM, we focus on 2D media, finding 
properties for the square, triangular and hexagonal WM 
(named respectively SWM, TWM and HWM in the fol- 
lowing). Such properties allow to calculate the bandwidth 
of a signal produced by a WM working at a given sam- 
pling frequency, once its geometry and the density of its 
junctions have been determined. 

The paper is structured as follows. Section || provides 
some background material on waveguide meshes (and their 
interpretation as finite difference schemes) and multidimen- 
sional sampling lattices. Section III illustrates the spatial 
sampling efficiency of the three waveguide mesh geome- 
tries for signals having circular spatial band shape. In 
Section IV we explain how the critical spatial sampling 
affects the choice of the temporal sampling frequency in 
non-aliasing conditions. In Section [v|, the computational 
performances of the three geometries are compared under 
critical sampling conditions. 

II. Background 

A. Digital Waveguides and Waveguide Meshes 

An ideal one-dimensional physical waveguide can be 
modeled, in discrete time, by means of a couple of par- 
allel delay lines where two wave signals, s+ and s_, travel 
in opposite directions. Such a structure, based on spatial 
sampling (with interval D) and time sampling (with in- 
terval T), is called a digital waveguide [Q, p3). If wave 
propagation in the physical medium is lossless and non- 
dispersive with speed c = D/T, no error is introduced by 
the discrete-time simulation as long as s+ and s_ are band 
limited to a band B = (2T) _1 . In this case, the signal 
s(x,t) along the physical waveguide can be reconstructed 
with no aliasing error from samples of the wave signals: 



s(mD,nT) — s+(mD,nT) + s—(mD,nT) 



(1) 



N digital waveguide terminations can be connected by 
means of a lossless scattering junction 6), || . Preservation 
of the total energy in the form of Kirchhoff's node equations 
leads to the scattering equation 



N 



N 



,N 



(2) 



k=l 



which allows for calculating the outgoing wave signal to 
the i-th waveguide branch from the N incoming wave sig- 
nals Si+, . . . , Sj\r-H under the assumption of equal wave 
impedance at the junction for all the waveguides. 
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Fig. 1 

The SWM (a), the TWM (b) and and the HWM (c). — ,/,\ and | are digital waveguides, • are lossless scattering junctions. In 
(c), seven lossless scattering junctions separated by two digital waveguide branches are marked with 0. 



A WM, as proposed by Van Duyne and Smith in 1993 || , 
[U, is obtained by connecting unit-length digital waveg- 
uide branches by means of lossless scattering junctions. 
For the simulation of uniform and isotropic multidimen- 
sional media, a few kinds of geometries, all corresponding 
to tiling the multidimensional space into regular elements, 
have been proposed: square ||, triangular Q], hexag- 
onal B for 2D media such as membranes; rectilinear |j) 
and tetrahedral 0, [|] for 3D media. The WMs which 
are considered in this paper (SWM, TWM and HWM) 
are depicted in figure [fl Among the proposed geometries, 
the HWM is peculiar because its 3-port lossless scatter- 
ing junctions exhibit two different orientations, and it can 
be interpreted as two interlaced TWMs (see the junctions 
marked with in figure [l].c). 

WMs introduce a basic relation, between signal 
s(Xj,nT), taken from a junction located at position Xj in 
space, and signals s(xj + Dfc, nT) , k = 1, . . . , N, taken 
from the N adjacent junctions connected to it, D meters 
far from position x.,-: 



s(x 3 ,nT+T) + s(x 3 ,nr- 



2 N 
k=l 



■D fc ,nT) (3) 



which is obtained from (Q). 

Equation (||) indicates that each WM behaves like a finite 
difference scheme the difference being that the former 
has a state lumped in the digital waveguides, and the lat- 
ter has a state lumped in the junctions. For the purpose of 
the analysis that follows, when the lossless scattering junc- 
tions have more than one orientation, as in the HWM, a 
difference equation such as (||) should be written as many 
times as there are orientations, each one using a proper set 
of vectors D l5 . . . , Djv @ ■ 

Following the lines of von Neumann stability analysis || , 
equation (|^) can be Fourier transformed with spatial vari- 
ables x and y, resulting in 

S(£x , £, y , nT + a g T) + S(£ x , £„ , nT - a g T) (4) 



where £ x and £ y are spatial frequencies, a g takes the value 
2 for the HWM and 1 for the other meshes, and bg is a 
geometric factor equal to 



b a 
h 



cos(2wD£ x ) + cos(2wD£ y ) 

2 2 / 

- cos (2irD£ x ) + - cos 27r£> 



2 tx + 2 Zy 



■ cos 2nD 



(5) 



cos 



■ cos 2tiD 



-— cos 2nD 
9 



2^ x 2^ 



2 



for the SWM, TWM and HWM, respectively. 

Solving equation (|J) as a finite difference equation in 
the discrete-time variable, the spatial phase shift affecting 
a traveling signal in one time sample is found to be 



&<Pg(€x>€y) = —— arctan 



4-6g 



(6) 



We can compare the propagation speed of a signal trav- 
eling along a WM, versus the propagation speed of a signal 
traveling along an ideal membrane. If we consider mem- 
branes where relation D — cT holds, meaning that signals, 
during a time period T, travel for a distance equal to the 
digital waveguide length, we can calculate the spatial phase 
shift of these signals, occurring during a time period: 



Aip = -2ttD£, , 



(7) 



where f = yCx + C y ■ By comparing (||) and (Q) , we find the 
ratio fcg between the propagation speed of a signal traveling 
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along a WM and along an ideal membrane 



• arctan 



(8) 



2wa g D£, u g 

Since A; is a non constant function of the spatial fre- 
quencies, WMs introduce a dispersion error, i.e., the signal 
traveling along a WM is affected by dispersion of its com- 
ponents. As a starting point, dispersion can be evaluated 
for spatial frequencies lower than the Nyquist limit£], that 
is, in the frequency domain: 



(9) 



This domain will be refined, according with the considera- 
tions to be presented in section III. 

Figure || shows plots of kg, kt and k^, where D has been 
set to unity. It can be noticed that the propagation speed 
decreases for increasing spatial frequencies. In particular, 
k s is maximum when £ x = £ y , and minimum for high val- 
ues of the spatial frequencies located along the main axes, 
suggesting that dispersion in the SWM does not affect the 
diagonal components traveling along it. On the contrary, 
the HWM exhibits the flattest dispersion error on the re- 
gion centered around dc. Finally, the TWM seems to have 
the most uniform behavior of the dispersion error. 

The propagation speed in all the WMs has a maximum 
at dc: 

fc g (0) = limfc g (0 = ^. (10) 

It is worth noticing that this value corresponds with the 
nominal propagation speed of a signal traveling along a 
finite difference scheme M. 

iFor the purpose of this paper, the Nyquist limit is denned as half 

B sample rate, both in time and space. In some previous works |7f , 
the simulations were considered valid up to a quarter of the sam- 
rate because in the square mesh the frequency response repeats 
itself after that limit. However, as it was pointed out in |7l, this is 
due to the fact that all transfer functions definable at any one junc- 
tion are functions of z~ 2 . This does not imply that the response to a 
signal having components up to half the sample rate will be aliased. 
However, for certain applications such as modal analysis of physical 
membranes, the frequency response at quarter of the sample rate is 
the definitive limit with the square mesh. 



B. Sampling Lattices 

The evaluation of the dispersion error does not give a 
complete description of the constraints holding when an 
ideal membrane is modeled using WMs. In particular, a 
method is needed for computing the signal bandwidth a 
WM is able to process. The theory of sampling lattices Q] , 
fLU , that we are briefly reviewing in this section, gives the 
background for characterizing WMs from this viewpoint. 

Let us sample a 2D continuous signal s over a domain L, 
subset of 1Z 2 , so defining a discrete signal s L (x) , xel. If 
L can be described by means of a nonsingular matrix L such 
that each element of the domain is a linear combination of 
the columns of L, the coefficients being signed integers: 



x = L 



»1 

Ho 



, mi € Z , U2 € Z , 



(11) 



then L is called a sampling lattice, and L is its basis. 
The number of samples per unit area is O] , p5[ : 



Vr = 



1 



dct(L) 



(12) 



as L contains information about the distance between ad- 
jacent samples. 

Sl, the Fourier transform of Sj,, is defined over K 2 and 
obtained by periodic imaging of S, Fourier transform of s. 
These Fourier images are centered around the elements of 
the lattice L* , which is described by the basis L~ T , inverse 
transposed of L. Notice that the denser the sampling of s 
is, the sparser the image centers are. In formulas, 



V 



1 



L * 



det(L- T ) 



det(L) 



1 



(13) 



Let us consider an ideal, unlimited membrane traveled 
by a spatially band limited signal s(x, y, t), and let us do a 
spatial sampling of the signal, at a given time. Whenever 
the spatial sampling defines a sampling lattice L, so that sl 
is defined, we can calculate Sl using the above results. The 
multidimensional sampling theorem fl4| ] — which, in brief, 
tells that if the Fourier images of s do not intersect one 
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with each other, s can be recovered from sl with no alias- 
ing error — indicates whether the chosen sampling scheme 
induces aliasing. 

Equation ( |l3| ) tells that for each choice of the sampling 
lattice geometry, the density T>l can be increased until the 
images do not intersect. Conversely, given s, there exists 
a sampling lattice capable to capture all the information 
needed to recover the original signal using the least density 
of samples. Clearly, it will exhibit the highest sampling 
efficiency. 

III. Sampling efficiency of the WMs 

Sampling efficiency will be calculated for signals having 
circular spatial band shape centered around the origin of 
the frequency axes^, with radius equal to B. Even if the 
analysis procedure does not depend on the shape of the 
spatial domain, the circular band seems to include all the 
signals occurring in practical applications. Indeed, the pro- 
cedure is independent of the system evolution — which in 
the WMs is controlled by equation (Q) — so that it applies 
to any model discretizing distributed systems, and where 
signal information can be located over a sampling lattice. 
This is the case, for example, in finite difference schemes. 

A. WMs and Sampling Lattices 

A SWM, having digital waveguides of length D s , corre- 
sponds to a sampling scheme over the lattice L S (D S ), de- 
scribed by the basis 



L S (A 



A 
A 



(14) 



A TWM, having digital waveguides of length £>t, corre- 
sponds to a sampling scheme over the lattice A (A), de- 
scribed by the basis 



L t (A) = 



A 




|_A 

#A 



(15) 



Notice that a triangular scheme is denser than a square 
scheme made with digital waveguides of the same length. 
This is confirmed by relation 



V 



L t (D) _ det(L s (A)) D 2 



V Ls{D) det(L t (£>)) 



f A 



2 

7i 



= ~m ■ ( 16 ) 



The description of an HWM having digital waveguides of 
length Ad again can be given in terms of sampling lattices, 
with some extra care. The HWM is obtained by subtract- 
ing a TWM, whose junctions lie on the sampling lattice 
At(Ai)j from a denser TWM having junctions on A(A)- 
Hence, the junctions of the HWM lie on 



A 



L h (D h )=L t (D h )\L T (D h ) 



The basis of A (A) ^ s 



L T (D h ) 



§A 

2 h 







(17) 



(18) 



y 



X 



D h 



X > 

\ 2 'V 
-* 1 — * — x 



Fig. 3 

The HWM, obtained by subtraction of TWMs. x are 

ELEMENTS BELONGING TO Lr^(D^), • ARE ELEMENTS BELONGING TO 

L t (D h ). 



Figure | shows the HWM over A(Ai)j obtained by sub- 
tracting the sampling lattice A (A) (whose elements are 
depicted with x) from A (A) (elements depicted with •). 

B. TWM vs SWM 

The image centers of the spectra belonging to signals 
traveling along a SWM and a TWM are respectively de- 
scribed by the basis matrices of A and A* : 



and 



L S (A 



L t (A 



l/A 






l/A 



^ o 



2 1 



(19) 



(20) 



Figure ^ shows Fourier images for the SWM (empty cir- 
cles in dashed line) and the TWM (filled circles), located 
around the origin of the frequency plane. It can be noticed 
that the square sampling scheme induces a square position- 
ing of the images and, consequently, a square tiling of the 
frequency plane, as emphasized by the square in dashed 
line. Similarly, sampling over a TWM results in a triangu- 
lar positioning of the images, thus producing an hexagonal 
tiling of the frequency plane, as shown by the hexagon lo- 
cated over the center. Such hexagonal tiling allows, as it 
appears quite evidently in the figure, to "pack" the images 
better than those coming from a square sampling. 

When both a SWM and a TWM critically sample the 
same signal (i.e. the filled circles touch each other without 
intersecting, and the same thing happens for the empty 
circles), it is easy to derive the following relation between 
the digital waveguide lengths: 



A 
A 



2 

7! 



1.134 



(21) 



2 This class of signals encompasses the signals obtained when a mem- 
brane is excited, in one or several points, by a single shot or by a 
sequence of shots, using a stick with an approximately round tip. 
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Fig. 4 

Domains of Fourier images in a SWM (empty circles in 

DASHED LINE) AND IN A TWM (FILLED CIRCLES), AND THE 
CORRESPONDENT TILING INDUCED BY THE RESPECTIVE GEOMETRIES 
(SQUARE IN DASHED LINE AND HEXAGON IN SOLID LINE). THE EMPTY 
AND THE FILLED CIRCLES HAVE THE SAME RADIUS AND TOUCH EACH 
OTHER WITHOUT INTERSECTING, MEANING THAT BOTH THE SWM 
AND THE TWM CRITICALLY SAMPLE THE SAME SIGNAL. 




Fig. 5 

Centers (x) of the Fourier images coming from a TWM 
defined by l^-d^)' and centers (•) of the fourier images 

coming from a sparser twm defined by 1/^(1)^). 
Subtraction of x from • gives the centers of the Fourier 
images (located on the filled circles) coming from an hwm 

defined by l^dfo). the respective tiling is grven by the 
triangles. both the hwm and the sparser twm critically 
sample the same signal. 



Under this condition we can relate the sample densities 
in the two geometries: 



V 



V 



det(L s (£> s )) 



V 



Ls(D s ) 



V 



Ls(Ds) det(L t (^D s )) 



f , (22, 



concluding that the TWM exhibits a better sampling effi- 
ciency relative to the SWM. In other words, a signal can be 
spatially sampled with a triangular geometry using 13.4% 
less samples per unit area. 

C. TWM vs HWM 

Said Sh and St the signals sampled over the junctions 
of the HWM and the TWM, respectively, and said sx the 
signal sampled over the lattice L^(Dif), we can define the 
zero-padded signals 



s T (x) 





s T (x) 



x € L h (D h ) 
x G L T (D h ) 



and 



s h (x) = s t (x)-s T (x) = 



Sh(x) 




x G Lh(Ai) 
x G L T (D h ) 



(23) 



(24) 



Since we cannot define a Fourier transform over Lj 1 (I?j 1 )> 
we consider S^, Fourier transform of (which is defined 
over Lt(-Dh))i as a description of s^ in the frequency do- 
main. S\y, according to its definition, can be obtained by 
subtracting St from St: 



(25) 



The result is shown in figure where some images of St 
are depicted (filled circles plus circle in dashed line). The 
elements of L^*(D]f) are marked with •, while the elements 
of Lt*(D]f) are marked with x. 

The HWM induces an hexagonal positioning of the 
Fourier images (corresponding to the filled circles). In fact, 
their centers are elements of a set which, again, can be 
defined by a subtraction between two sampling lattices, 
Lt*(.Dh) an d -^T*(-Dh)i which are reciprocal of Lt(D^) and 
Irx(-Dh); respectively. The consequent tiling geometry is 
triangular, as emphasized by the triangles depicted in fig- 
ure ||. 

It can be observed that image intersection does not occur 
in Lh(Dif) if and only if s is sampled without aliasing over 
Lx(-Dh)- I n other words, when 5*t exhibits superposition 
of its images, the same thing happens for the images of S^, 
and vice versa. 

Moreover, it must be noticed that equation 



L T (D h ) = L t (V3D h 



(26) 



holds between Lt and Lt, if a rotation is neglected. This 
relation, together with the considerations about image su- 
perposition made just above, allows us to say that an HWM 
does not perform better than a TWM, whose digital waveg- 
uides are y/3 times longer. 

Another interesting consideration comes out by noticing 
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that, since 



V L t (D h ) = det(L T (J> h )) 
■D LT( D h ) det(Lt(Dh)) 



= 3 



(27) 



then the TWM defined over the lattice ^(D^) is 3 times 
as sparse as the TWM defined over Lt(D^), thus (directly 
from the definition of L n ) 2 times as sparse as the HWM. 
Hence, the number of lossless scattering junctions per unit 
area of the HWM is twice as high as that of the TWM 
defined over Lx(-Dh)' with no benefits on the accuracy of 
sampling. This relates with the fact that the hexagons 
tiling the frequency plane (as the inner hexagon depicted 
in figure contain twice the information needed to recover 
s correctly: the 6 partial images included in the slices inside 
the central hexagon can be composed into 2 Fourier images. 

IV. Signal time evolution 

In this section, we discuss the consequences of critical 
spatial sampling on the temporal sampling frequency. 

First, let us consider the sampled version of a signal trav- 
eling at speed c along an ideal membrane that has been ex- 
cited by a bandlimited signal f(x,y) having Fourier trans- 
form equal to F(^ x , £ y ). It can be shown (see Appendix A) 
that a time sampling frequency 



2c max {£} 



(28) 



is required to recover the original signal correctly. 

When such a signal has a spatial circular band shape, 
thus belonging to the class seen in section [II, the relation 



max j£) = B 
t-\F{t)\ita 



holds, so that equation 



simplifies into 



F, = 2cB . 



(29) 



(30) 



Relation (|30|) has an immediate physical interpretation: 
waves having wavelength equal to l/B, propagating at 
speed c along a medium, exhibit a temporal frequency equal 
to cB. In order to preserve information in their sampled 
versions, they must be sampled above twice their frequency. 

The value given by equation ( |30| ) , if used as the sampling 
frequency of a 2D-resonator model (realized by means of 
WMs having critical waveguide lengths), causes inaccurate 
positioning of the modal frequencies. Assumed that disper- 
sion (see section II- A) cannot in general be eliminated, as 
it is a consequence of the finite number of directions a 2D 
signal can propagate along a WM0, the propagation speed 
can at least be set to its physical value in low frequency. 
This can be done by simply rescaling F s according to the 
geometry. 

Recalling the procedure leading to equation (0), we can 
reformulate the spatial phase shift of a signal traveling 

3 Savioja and Valimaki jn| have pointed out that a signal, coming 
out from a TWM, can be frequency warped to reduce the disper- 
sion error. This is made possible by the fact that k^(^ x ,^ y ) can be 
approximated with a single- variable function fct(£). 



along the ideal membrane, during a time period T = l/F s . 
Hence, we find the following expression for the ratio (|§|), 
reformulated using for each geometry its respective critical 
waveguide length: 



1 



2ira g D£, 



arctan 



(31) 



where b g has the structure of 6 g (see equation (^)), but 
the waveguide length D has been replaced by its critical 
counterpart Dg. 

Recalculation of the limit ([h]) gives: 



Mo) = ^ g (o) = ^§, 



(32) 



that is, considering D s as the reference waveguide length, 



fcs(0) 



1 

71 



fct(0) 



V2 
V3 



fch(0) 



V2 
3 



(33) 



This result shows that signals, under conditions of crit- 
ical sampling, propagate in the WMs at different speeds, 
according to the geometry. The dispersion ratios can be set 
to unity at dc by using the temporal sampling frequencies 



Fs, g = ~ 



1 



fcg(0) 



-F s 



(34) 



For the purpose of sec. the dispersion ratio of a critically 
sampled 2D medium, adjusted to be one at dc, is called 

V. Performance 

In order to compare the three geometries under critical 
sampling conditions, we show in figure || the contour plots 
of kg, for a nominal spatial bandwidth B — 1/2 (corre- 
sponding to the circle in the figure). 

It can be noted that the behaviors of the SWM and the 
TWM are not dramatically different in terms of average dis- 
persion error: dispersion stands quite below 20% for most 
spatial frequencies in both the geometries. However, the 
TWM exhibits a more uniform behavior, and this unifor- 
mity can be exploited using frequency warping JTo[ | . Con- 
versely, dispersion in the HWM stays below 10% almost ev- 
erywhere, thus indicating that this geometry has the most 
uniform propagation speed under these test conditions. 

The computational cost and the memory requirement in 
the different geometries can be calculated under the same 
conditions. They are based on equation (^) which shows 
that each iV-port lossless scattering junction requires 27V 
operations to compute the wave signals coming out from 
the junctions, to be stored into N locations belonging to 
the adjacent digital waveguides. These operations amount 
to 2N — 1 additions, and 1 multiplication (which can be 
replaced by a bit shift in fixed-point implementations of 
the SWM). 

Table | summarizes the performance, when the reference 
sampling rate of the square mesh F ss has been set to a 




TABLE I 

Performance of the geometries in terms of memory requirement and computational cost (multiplications are pure bit-shifts 
in the SWM). Two implementations are considered: as a waveguide mesh and as a finite difference scheme. 





Waveguide Mesh 


Finite Difference 
Scheme 




Sq. 


Tr. 


Hex. 


Sq. 


Tr. 


Hex. 


additions per junction 


7 


11 


5 


4 


6 


3 


multiplications per junction 


1 


1 


1 


1 


1 


1 


memory locations per junction 


4 


6 


3 


2 


2 


2 


density of junctions 


1 


0.866 


1.732 


1 


0.866 


1.732 


density of memory locations 


4 


5.2 


5.2 


2 


1.732 


3.464 


sample rate 


1 


0.866 


1.5 


1 


0.866 


1.5 


additions per unit time and space 


7 


8.25 


12.990 


4 


4.5 


7.794 


multiplications per unit time and space 


1 


0.75 


2.598 


1 


0.75 


2.598 



nominal value equal to unity. Numbers are given for two 
implementations: as a waveguide mesh (with memory in 
the waveguide branches), and as a finite difference scheme 
(with memory in the junctions). 

The numbers of operations and memory locations per 
junction for the WM follow directly from equation (||). 
The numbers of operations and memory locations per junc- 
tion for the finite difference scheme follow directly from 
equation (|j). The densities of junctions result from equa- 
tions ([H]) and ( p6| ) . The densities of locations are obtained 
by multiplying the third row times the fourth row. The 
sample rates are a consequence of equation (|i|) . 

The results of table || show that the triangular finite dif- 



ference scheme uses the least quantity of memory. Among 
WM implementations, the SWM is the most efficient in 
terms of memory occupation. 

Finally, the number of additions (multiplications) per 
unit time and space results by multiplication of the num- 
ber of additions (multiplications) per junction, density of 
junctions, and sample rate. Once again, the SWM has the 
least density of operations when a WM model is required. 
Conversely, the square and the triangular geometries have 
about the same computational requirements in a finite dif- 
ference implementation. 
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A. Numerical example 

Let us design a WM, capable of modeling in real-time 
an ideal 2D round resonator, of radius r = 0.1 m, where 
waves propagate at a speed c equal to 130 m/s. Let the 
signal contain information up to a frequency / = 10 kHz. 

The spatial bandwidth B of the signal traveling along 
the resonator is: 



/ 10000 Hz _ t 

B = — = — — — = 76.923 m 

c 130 m/s 



(35) 



and consequently the critical waveguide lengths required to 
model the signal, in the respective geometries, are: 



A 



1 

2B 
1 



V3B 
1 

3B ~ 



6.5 mm ; 
= 7.5 mm 
4.3 mm . 



(36) 



This means that the numbers of junctions — iV s , iVt and 
Nfr, respectively — needed in the three models areFI 



Furthermore, when processing signals of the same band- 
width, the triangular waveguide mesh does not exhibit a 
computational and memory load much larger than the most 
computationally-efficient waveguide mesh — the square 
mesh — such that it can be considered, for the unifor- 
mity of its dispersion error, a good choice both in terms of 
simulation errors and computational cost. 

The analysis presented in this article may be extended 
to 3D media, comparing the three geometries which most 
directly correspond to the waveguide meshes here reviewed: 
3D rectilinear, dodecahedral, and tethraedral. 
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Appendix A: Sampling of a signal traveling 

ALONG AN IDEAL MEMBRANE 

An ideal membrane establishes a relation between the 
spatial Fourier transforms of the signal s traveling on it, 
taken in correspondence of two times, t and t: 



N s 
N t 
N h 



744 



2 irr 2 



3 V3 D 2 
2 u h. 



= 645 



= 1308 



(37) 



The time sampling frequencies needed to have kg(0) = 1 
are, from equation 



Fs,s 
F s ,h 



= 2V2f = 28.285 kHz 



(38) 



-f = 24.495 kHz ; 

2-4=/ = 42.427 kHz . 
V2 

VI. Conclusion 

A novel wave propagation model has recently been in- 
troduced for the simulation of isotropic multidimensional 
media. It makes use of structures called waveguide meshes. 
Even if most of the waveguide mesh properties have already 
been understood, there was lack of literature about their 
performance from a signal-sampling viewpoint. 

In this paper, some properties of the most common 2D 
waveguide meshes — square, triangular and hexagonal — 
related with the bandwidth of the signal traveling on them, 
have been inspected. In particular, it has been shown that 
the triangular waveguide mesh is capable of processing a 
larger bandwidth than the square or hexagonal waveguide 
meshes having the same digital waveguide lengths. 

4 The values result by calculating the number of small areas, each 
one being associated with its own lossless scattering junction, which 
tcsscllate the resonator. 



(39) 



This relation has the following interpretation: each spatial 
component of the signal, traveling along a distance equal 
to c(t — t) during time [t, i), has no magnitude variation, 
and has phase variation equal to — c(t — £)£ 

Let us excite, at time t, the ideal membrane with a spa- 
tially band limited signal f(x, y) having Fourier transform 
F(£ x ,£ y ). The application of ( |39| ) gives the spatial Fourier 
transform of the signal on the membrane: 







t < t 



F(C x ,tiy) , t>t 



r (40) 



From this relation, we can calculate the critical tempo- 
ral sampling frequency required to sample, without loss of 
information, a signal traveling on an ideal membrane. In 
fact, the Fourier transform of (fhj), 



/+oo 
-OO 



(41) 



has a magnitude equal to 

\S&c,ty,f)\ 



+ 00 



+ 00 



+oo 



J2-Kc{t-t)i 



F(£ x ,ty)e- j2 * ft dt 



(42) 



\F&,Zv)\ 



e - j2 nt{f-ci} dt 
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In particular, when t 

\S(&,Z v ,f)\ 

= \F(C x ,Q\Jim 



\F(Z*,Zv)\ 



t— y — oo 



+ OC 



e -,-2*t{/-ce> dt 



(43) 



This means that some spectral power exists for any tempo- 
ral frequency (equal to c£) associated to a spatial compo- 
nent of frequencies (£ x , £ y ) excited in the membrane. Invok- 
ing the sampling theorem the critical time sampling 
frequency F s results to be equal to 
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[3] 
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